Active optical frequency standard using sequential coupling of atomic ensembles 
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Recently, several theoretical proposals adressed the generation of an active optical frequency stan- 
dard based on atomic ensembles trapped in an optical lattice potential inside an optical resonator. 
Using atoms with a narrow linewidth transition and population inversion together with a "bad" 
cavity allows to the realize the superradiant photon emission regime. These schemes reduce the 
influence of mechanical or thermal vibrations of the cavity mirrors on the emitted optical frequency, 
overcoming current limitation in passive optical standards. The coherence time of the emitted 
light is ultimately limited by the lifetime of the atoms in the optical lattice potential. Therefore 
these schemes would produce one light pulse per atomic ensemble without a phase relation between 
pulses. Here we study how phase coherence between pulses can be maintained by using several 
inverted atomic ensenbles, introduced into the cavity sequentially by means of a transport mecha- 
nism. We simulate the light emission process using the Heisenberg-Langevin approach and study 
the frequency noise of the intracavity field. 

PACS numbers: 06.20.fb; 42.60.Jf: 42.60.Rn 
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I. INTRODUCTION 

The concept of an active frequency standard in the op- 
tical domain was recently proposed and studied by sev- 
eral authors [ll-ll]- The main idea is to create a "super- 
radiant" laser operating deep in the bad cavity regime 
where the decay rate k/2 of the cavity field significantly 
exceeds the linewidth jab of the lasing transition. Then 
the fundamental linewidth D of the steady-state operat- 
ing laser can be described by the generalized Schawlow- 
Townes formula that can be written as 
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(1) 



in the practically interesting range A ^ {'^ab + 0- 
[T]| . Here A = Wa — Wc is the difference between the 
cavity resonance frequency loc and the frequency Wa of 
the lasing transition, g is the coupling coefficient, Mao is 
the steady-state occupation of the upper lasing level. The 
frequency uj of the emitted radiation is connected with 
the resonance frequency lOc of the cavity and with the 
frequency Ua of the atomic transition via cavity pulling, 
namely 
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In the bad cavity regime, any fluctuation (thermal or 
mechanical) of the cavity length has much less influence 
on the spectrum of the output light than in conventional 
lasers, operating in the good cavity regime. Fluctuations 
of the cavity length are currently the main limiting factor 
for the performance and even more for the transportabil- 
ity of modern passive optical frequency standards. An 
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active optical frequency standard could help to overcome 
these limitations. 

Two main approaches towards an active optical fre- 
quency standard are proposed today: the optical lattice 
laser [ll-Q and the atomic beam laser The first ap- 

proach suggests to use trapped atoms with narrow optical 
transitions (such as -'^Sq -O- '^Pi, {i = 1,2,3) transitions 
in divalent atoms) confined to the Lamb-Dicke regime in- 
side an optical lattice potential as a gain medium to build 
the laser. The necessary population inversion can be pro- 
vided by additional repumping fields, coupling the lower 
lasing state with some higher levels from which the atoms 
decay to the upper lasing level (3- level laser scheme). The 
general operation principle of such a optical lattice bad 
cavity laser has been demonstrated experimentally us- 
ing a Raman system to mimic a narrow linewidth optical 
transition [l2j . 

The main limitation to the performance of this optical 
lattice laser is the limited lifetime of the atoms in the op- 
tical trap. A laser whose gain medium is a single ensem- 
ble of trapped atoms can not keep the phase longer than 
the atom trap lifetime. Another problem emerges from 
the pumping laser which can shift the frequency of the 
lasing transition. Although the authors of [3 claim that 
ac Stark shift-induced fluctuations of the atomic transi- 
tion frequency are negligible at the mHz level for repump- 
ing rates of the order of 10"^ s~^ and for a pumping laser 
intensity-stabilized to 1 %, also the frequency of the laser 
depopulating the low lasing level needs to be quite well 
stabilized. In Q, it is proposed to use a four- level lasing 
scheme with a non-stable lower lasing state to avoid the 
pumping- induced shift, however this scheme still suffers 
from the problem of limited trap lifetime. 

The atomic beam approach allows to overcome the 
atom lifetime limitation. Here, a continuous beam of ac- 
tive atoms, previously pumped to the upper lasing state, 
passes through the cavity. The theory of laser generation 
for this system was presented in Q and the generalisa- 
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tion to Ramsey-type interaction between the atoms and 
the intracavity hght was performed in In the au- 
thor claim that laser radiation with 0.5 Hz linewidth and 
120 nW power can be attained with a beam of hot ®®Sr 
atoms prepmnped to the "^Pi upper lasing state (natu- 
ral linewidth 7 kHz), with a cavity of 800 /^m waist and 
k/2 = 11 MHz. The flux required for the lasing is 10^^ 
atoms per second. At the same time, the theory de- 
veloped in [1] does not take into account the motion of 
the atoms along the cavity axis while passing through 
the resonator. This motion will lead to additional line 
broadening and a first-order Doppler shift. The latter 
will fluctuate together with the mean transversal atomic 
velocity which depends on the environment, thermal fluc- 
tuation of the nozzle etc. Accurate compensation of all 
these effects appears to be a very challenging task. In Q , 
a similar idea is proposed for the creation of an extreme 
ultraviolet laser using mctasablc noble gas beams. 

We note that for lasing transitions with longer wave- 
lengths (such as infrared transitions in some molecules) 
it may be possible to combine the Ramsey-type lasing 
scheme Q with mechanical blocking of the particles that 
would pass through the interaction regions with different 
phases of the field, similarly to [T3| . 

In the present work we investigate an intermediate, al- 
ternative approach to the schemes presented above, con- 
sisting in dynamically moving trapped atomic ensembles 
into (and out of) the cavity. These ensembles will be pre- 
pared in the upper lasing state outside the cavity, circum- 
venting perturbations due to ac Stark shifts. Transport 
will be realized by means of a moving red-detuned ID 
optical lattice, transport of more than 2 ■ 10^ ultracold 
atoms over a distances of up to 20 cm with fast trans- 
port velocities of up to 6m/s and strong accelerations 
of up to 2.6 • lO^m/s'^ has already been demonstrated 
experimentally [l3|- To confine the atoms to the Lamb- 
Dicke regime inside the cavity, an auxilary blue-detuned 
stationary optical lattice is used in our scheme. Both, 
the red-detuned and blue-detuned optical lattice should 
have "magic frequencies" [l^ providing equal light shift 
for both states of the lasing transitions. 

The operational sequence for two ensembles consists 
of 4 stages. In the first stage, one atomic ensemble is 
positioned inside the cavity and starts to emit light into 
the cavity mode whereas a second ensemble is trapped 
outside the cavity, pumped into the upper lasing state 
and transporting towards the cavity, see Fig. [TJ In the 
second stage, when a significant part (but not all) of the 
atoms in the first ensemble are transferred to the ground 
state, the second ensemble is introduced into the cavity. 
The atoms of the second ensemble start to emit photons 
into the cavity mode mainly via stimulation emission, 
maintaing the optical phase. In the third stage the first 
ensemble is extracted from the cavity and a new inverted 
ensemble is prepared (either by loading new atoms or by 
repumping) while the second ensemble emits the light. 
The fourth stage is identical to the second one up to per- 
mutation of the ensembles. To avoid the time-dependent 
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FIG. 1: Schematic of two atomic ensembles (black dots) in 
two moving optical lattices during the first stage of the las- 
ing process. The first atomic ensemble emits light into the 
resonator mode while being slowly pulled through the cavity 
waist, the second one is prepared and transported towards 
the cavity. 



second-order Doppler effect, one should keep a constant 
velocity of the moving optical lattice while atoms cross 
the cavity waist. This process can be repeated as many 
times as necessary to keep the phase of the laser field. 

In the following, we model the lasing process using 
the Heisenberg-Langevin approach and study the phase 
stability of the emitted light. 



II. MODEL 

A. Operator Heisenberg-Langevin equations 

We employ a model based on the set of quantum 
Langevin equations following 0, [l^ . One or two ensem- 
bles Ai and A2 of homogeneously broadened two- level 
atoms with transition frequency uJab are held in a cavity 
characterized by a cavity damping constant k. The fi- 
nesse of the cavity is assumed to be sufficiently high to 
justify the mean field approximation. The atoms interact 
with the radiation field of a single excited mode of the 
cavity which we consider as a plane wave with frequency 
LOc- By adding up the individual atomic operators, we 
define the macroscopic atomic operators 

N^it) = (3) 
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where x = 1, 2 is the index of the corresponding atomic 
ensemble, al = (|6)(a|p, = i\a}{a\y, - 
a) and |6) are the upper and lower lasing states respec- 
tively. The equation for the photon annihilation operator 
a is: 

a = -'la + g (Af^ri(i) + Af2r2(t)) + F^{t), (4) 

where r^{t) describes the time dependence of the cou- 
pling of the xth atomic ensemble with the cavity mode, 
F-y is the Langevin force. Here we suppose for the sake 
of simplicity that all the atoms have the same coupling 
coefficient g and that LOab ~ '^c- The last approximation 
is valid if \ujab - Wc| < -|- 7ab). 

In our model we suppose that when a certain atomic 
ensemble starts to interact with the cavity field, all the 
atoms of this ensemble have been pumped into the state 
a). During the interaction time, the set of equations for 
the macroscopic atomic operators ([3]) of xth ensemble is: 

h = laN^ + 9^x{t) + MX+a) + F^{t). (5) 

where 7a is the spontaneous decay rate of the upper lasing 
state, 7ab is the decay rate of atomic coherence, and Fm, 
Fa and f}, are the Langevin forces. Here we neglect loss 
of atoms from the trap during the interaction time. It 
is easy to show that Fa = —Fi, in this case. The mean 
values of the Langevin forces are @ : 

(P^it)) = (P.^t)) = {P>^,{t)) = (P.it)) 0, (6) 

Correlation functions of the Langevin forces correspond- 
ing to the intracavity field are 

(^P+it)P^{t'))^KnT5{t~t'), (7) 

= + mt - t'), (8) 

{P+it)P+it')) = {P,it)P,it')) = 0, (9) 

where fix is the number of thermal photons in the cavity 
mode. We suppose that the temperature is low enough 
to set fiT = 0. In this case the correlations between the 
Langevin forces for a specific atomic ensemble are 

{P^it)P^it')) = 7a (Nm) Sit - t'), (10) 

(F>^+(t)F^,(t')) = (27a. - 7a) (N^it)) Sit t'), (11) 

(p^it)Pi+it')) = (27a. (N^it)) + 7a (N^it)) ) 

xSit-t'), (12) 



{P^+it)P^+it')) = [Plit)Fl,it')) = 0, (13) 

l^PmFltit')) - {PlHnP^it)) = 0, (14) 

{Piit)Pitit')) = laiM^^itmt - t'), (15) 

{P]tAt')P^it))^7a{MHt))Sit-t'). (16) 

Langevin forces corresponding to different atomic ensem- 
bles are uncorrelated. 



B. Equivalent c-number stochastic Langevin 
equations 

As a next step we introduce the stochastic c-number 
Langevin equations which are equivalent to the oper- 
ator Langevin equations dH) and ([S]). To estabilish a 
unique relation between operators and c-number vari- 
ables we have to define the correspondence between prod- 
ucts of c-numbcrs and products of operators. In [16| and 
in a number of subsequent papers 0, |^, fll| . the nor- 
mal ordering a+, M+, Na, Nb, M, a is chosen and the 
c-number Langevin equations were derived in such a way 
that the equations for first and second moments of c- 
number variables A, A4, Ma, Afb are identical to corre- 
sponding equations for normally ordered operator vari- 
ables a, M, Na, Nb- This choice leads to a redefini- 
tion of correlation functions for c-number Langevin forces 
J- a, J'b, J'M, J'M*- In our case this choice of normal or- 
dering leads to self-contradicting correlation functions, 
see Appendix A for details. Moreover, the choice of nor- 
mal ordering for atomic operators M, Na, M+ is ambigu- 
ous, in contrast to the choice of normal ordering of field 
operators a and a^. 

To overcome these difficulties we use here the follow- 
ing rules: 1) The equations for first moments of operators 
and c-number values should be identical; 2) the equations 
for second moments of c-number field variables A and A* 
should be identical to the corresponding mean values of 
normally ordered field operators d and a+; 3) the equa- 
tions for second moments of c-number atomic variables 
Mx, A/"^, A^* should be identical to corresponding equa- 
tions for mean values of symmetrized products of operator 
variables M^, N^,M^^. Using these rules one can find 
the correlation functions for c-number Langevin forces in 
a similar way as in [T6| : 

(mt)) = {:F^{t)) = {:FUt)) = (j-^(t)) = 0, 
{T;{t)T,it')) = {T;{t)T;{t')) (i?) 

= {J'-,it)J^^{t')) = 0, 

and 

{^^it)^f{t'))=2Vl,Sit-t')Sx^, (18) 

where the diffusion coefficients 2?^; are 

2Pl = 7a(A/?(t)), (19) 
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= ^I^M'M' = 0, (20) 
^I^MM'^^abM^, (21) 

225^A4 = Y-^x, (22) 

22??A,. = Y^l- (23) 

Here = A^x +A/;>^. 

Equations Q and ^ are already written in proper 
order, so that it is easy to get the equations for the cor- 
responding c-number variables (for the sake of brevity, 
the argument t is omitted here): 

A= -'^A + giMiri+M2T2), (24) 



= -jabM^ + gT^ {2N^ -N^)A + IF^, (25) 
= -laN^ - gV^ {A*M^ + AM*^) + F^. (26) 



In the bad-cavity regime, the field variable A follows 
the atomic variables adiabatically: 



A = ^{MiT^+M2T2)- (27) 
It yields the full set of equations for the atomic variables 



--iabMl{t)+G 



M2{t) 

Ml{t) = -^Mt)-G 



Mi{t) -Tlit) + M2{t) ■T,{t)T2{t) ■ [Mlit) 



Mi{t)Ml{t)-Ti{t) + 



Mi{t)M*2{t) + Ml{t)M2{t) 



ri(i)r2(i) 



labM2{t) + G M2{t) ■Vl + Mi{t) ■T^{t)T2{t) ■ [Af^{t) 



M2{t)M*2{t)-Tl{t) + 



M2{t) 



Mi{t)M*2{t) + Ml{t)M2{t) 



T^{t)T2(t) 



Kit), 



Kit), 



(28) 
(29) 
(30) 
(31) 



r 



where G = Ag"^ / n. Using equations ([25)) — (PT|) and 
correlation functions for Langevin forces — p3)) . we 
simulate the lasing process. 



C. Simulation of c-number Langevin forces 



In the present work we solve the set of equations ([28]) - 
(PT|) using the simplest Euler algorithm. Namely, intro- 
ducing vectors x = |^A^i, 7V^, A^2, rewrite 
HI) as 



x{t)=r{t,x) + f{t), 



(32) 



where r(t,x) is the deterministic ("regular") part of the 
right-hand side of equations ([25]) - (PU)) . and f{t) = 
(Tl^{t),T^{t),J'li{t),J^^{t)) is a stochastic part. Then 
the algorithm is 



x{ti+i) = x{U) + r{U,x) + f{U) 
where At = U+i — ti, and 



t,+l 



1 



f{t)dt 



(33) 



(34) 



is a set of random variables with zero mean and with 
correlations 



fk{ti)fm{t]) 



(35) 



Here Vkm is one of the diffusion coefficients - (|23p . 
Note that stochastic terms corresponding to different 
atomic ensembles arc independent. 

It is convinicnt to introduce two real column vectors 
¥^{U) and ¥^{t,) such that 



r^it,) 



( mt) \ 

Rc[J-^,(t.)] 

vim[^^,(t,)] ; 



(36) 



Then the correlators psp can be rewritten in matrix form 
as 



{¥^U)mt,))'^) 



At 



(37) 



where the covariance matrix D-*^ can be derived from H 
^ (I35|): 



1 



\iHMV^ 



(38) 



(Here and below in this section an index x ^-nd an argu- 
ment ti are omitted for the sake of brevity) . The matrix 
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D is positive semi-definite according the Sylvester's cri- 
terion. Indeed, the determinant of this matrix is 



4At 



MMalab-MM*^ 



> 0. 



because \M\ < MaMb < MaN, and 7a& > 7a/2. 
Consider a column vector 

A = -F, 



(39) 



(40) 



where X = {X''^^\ X^'^\ X^'^")) is a matrix constructed of 
the normalized eigenvectors of ED corresponding to the 
eigenvalues Ai, A2, A3. Correlations of elements of A can 
be found as 



{Ak{ti)A„i{tj)) = Si-jSkmK. 



(41) 



Therefore the elements of A are independent random val- 
ues with zero means and given dispersions. These values 
can be easily simulated, and the Langevin forces can then 
be calculated by the inversion of PO)) . 



III. SIMULATION 

In this section we specify our model and present the 
results of simulations of the lasing process. Consider the 
functions Ti{t) and ^2{t) describing the time dependence 
of the coupling between an atomic ensemble and the cav- 
ity mode. We suppose that this coupling is provided by 
the motion of the optical lattice confining the atoms. It 
seems to be necessary to keep a constant velocity of the 
atomic ensemble during its motion through the cavity 
to avoid the time-dependent 2nd order Doppler effect. 
Usually the transversal distribution of the electric field 
is Gaussian. For this reason, we take the functions Ti{t) 
and r2 {t) of the form 



r^(<) = exp 



18(t-t?)2 



y2 



(42) 



in ith lasing cycle. Here tf is the time when the xth 
ensemble crosses the symmetry axis of the cavity, and T 
is the time for moving the atomic ensemble through the 
cavity between two points where the amplitude of the 
intracavity field is e^^^^ ~ 1% of its value (on the cav- 
ity axis). We neglect any coupling outside these points. 
Therefore T = ti + 2t2, where ti is the duration of stages 
1 and 3, and n is the duration of stages 2 and 4, see 
Fig. [2 

We now simulate the lasing process and study the de- 
pendence of intracavity field phase fluctuations on the 
temporal parameters of the lasing cycle, summerized in 
the ratio 



R 



(43) 



which describes the temporal overlap of different atomic 
ensembles inside the cavity. It follows from general con- 
siderations, that the lower R (more overlap), the better 




FIG. 2: Time-dependence of atom-cavity coupling expressed 
by the coupling functions ri(t) and T2{t) for different stages 
of the lasing cycle. 



the polarization of the "old" ensemble will be inherited 
by the "fresh" one. An ideal system would be a contin- 
uous delivery of inverted atoms trapped in the optical 
conveyor belt, however such a device remains to be re- 
alized experimentally. If the experimental setup traps 
the atoms in a pulsed regime and if it contains only two 
moving optical lattices, then the atomic ensemble should 
be loaded into the optical lattice, prepaired in the upper 
lasing state and transferred into the cavity during the 
time Tl . Therefore ri (and R) can not be made arbitrary 
small. This difficulty could in principle be overcome by 
scaling to more than two moving optical lattices in one 
setup, however significantly enhancing the experimental 
complexity. Therefore it is interesting to find the opti- 
mal range of R where the phase transfer between two 
ensembles still works. 

In our simulation we use the following parameters for 
the atomic ensemble and the laser cavity: k = 10^ 
s'^ (this value can be attained in a cavity with a fi- 
nesse of about 5-10'' and a cavity length about 10 cm), 
7a = lab ~ 5 • lO^'^ s^^ (this value is typical for the spon- 
taneous decay rate of the state of divalent atoms, see 
[13, [3), 9 = 2.5 s~^ (this value corresponds to a cavity 
mode waist of about 1 mm, the wavelength of the lasing 
transition is about 700 nm and the cavity length is about 
10 cm), and an average number of atoms in one ensemble 
J7 ~ 2-10^. An important feature of any realistic prepa- 
ration and trapping of the atomic ensembles is that the 
number of trapped atoms fluctuates. We suppose that 
the dispersion of the number of atoms in one ensemble is 
10%. 

For all lasing stages we simulate the lasing process ac- 
cording the method described in Section |lT] using the 
functions (|42|) . At the start of stage 2 and stage 4 of 
the lasing process, we chose the number Af of atoms in 
the second and the first ensemble respectively, including 
random fluctuations. We suppose that atoms enter the 
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FIG. 3: (Color online) Phase diffusion coefficient versus the 
duration of the lasing cycle for different values of the ensemble 
overlap parameter R. 



cavity mode in the upper lasing state |a) with zero polar- 
ization. We then trace the intra-cavity field variable A 
over 1500 lasing cycles. This process is repeated several 
times which allows us to study the Allan dispersion cr^ 
of the phase </) of the laser field: 



1 



(44) 



where and are the phases averaged over ith and 
(i + l)th adjanccd time intervals, each of which has a du- 
ration of At. In our case o'^(At) cx At. as it is typical for 
a random phase walk. This process can be approximately 
described by the equation 



dt 



fit), (/(t))=0, {f{t)f{t'))^2DS{t-t'). (45) 



It is easy to show that the phase diffusion coefficient 
D is connected with the Allan dispersion as 



aliAt) ^ —At 



(46) 



So one can easily find the phase diffusion coefficient D 
usind the calculated Allan dispersion. 

Dependences of the phase diffusion coefficient D on the 
duration 2(ti +T2) of the lasing cycle for different values 
of R are presented in Fig. [3l To discuss these results, let 
us consider the time dependences of amplitude |^| and 
phase (j) of the intracavity field A for some particular re- 
alizations, sec Fig. m We see that if the lasing cycle is 
too short, the atoms have not enough time to develop a 
significant polarization. Therefore almost all the atoms 
leave the cavity before they can transfer their energy to 
the cavity mode, and the intracavity field is extremely 
weak, see Fig. |4] (a). Weakness of the intracavity field 
leads to a comparatively large contribution of Langcvin 
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FIG. 4: Behaviour of amplitude and phase of the intracavity 
field in a particular realisation over several lasing cycles for 
R = 0.25 andri-|-T2 = 0.1s (a), 0.25 s (b), 0.6s (c). Here non- 
zero to = 12 s was taken to show this behaviour in established 
regime 



forces which in turn leads to a large phase diffusion co- 
efficient. On the other hand, if the lasing cycle is too 
long, the atoms of the first ensemble transfer all their en- 
ergy to the intracavity field during the stage 1, see Fig. 3] 
(c). Therefore at the beginning of the stage 2, when the 
"fresh" ensemble is introduced into the cavity, the first 
ensemble has already lost the polarization. Of course, 
any transfer of polarization and passing on of the phase 
is impossible in this case. For the optimal duration of the 
lasing cycle, atoms have enough time to develop signifi- 
cant polarization, and this polarisation survives until the 
introduction of the next atomic ensemble into the cavity, 
see Fig. 111(b). 

We obtain optimal values for the lasing cycle of about 
0.4 - 0.6 s length for the experimental parameters intro- 
duced above. This result can be understood intuitively 
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FIG. 5: Power spectral density of the intracavity field for 
n = 0.0625 s, T2 = 0.1775 s 

by comparing it to the characteristic duration tp of a 
single superradiance pulse of an atomic ensemble cou- 
pled to the cavity, see Appendix B. For a single ensemble 
(TV = 2 • 10^ G = 2.5 • 10""*) we estimate tp ~ 0.1s. In 
the dynamic implementation discussed here, this time 
should be increased following the time dependence of 
r(i). Because one lasing cycle consists of two superra- 
diation pulses, one understands that the optimal tim- 
ing of these pulses corresponds to about 2 - itp. Also 
we note that the optimal lasing cycle duration decreases 
with smaller R because the simultaneous coupling of 
two atomic ensembles with the cavity mode results in 
a shorter radiation pulse. 

Finally, let us consider the one-sided power spectral 
density 

Sa{v) = ^ — 

oo 

= A j (Re(^(0)^*(T)))cos(27rz/)dT, (47) 



where z/q = 27raJab- The power spectral density for 
Ti + T2 = 0.25 s (which corresponds to a total duration of 
one lasing cycle of 0.5 s), R = 0.25, averaged over 60 re- 
alizations, each with 1500 cycles, is represented in Fig.[S] 
One can observe a high central peak corresponding to 
the main lasing component with the frequency vq, and a 
number of sufficiently lower sideband peaks at frequencies 
corresponding to multiples of the repetition frequency of 
the transport cycles (2 Hz in our case). 



dard by sequental coupling and decoupling of prepumped 
atomic ensembles with the optical cavity mode. We stud- 
ied the phase diffusion of the intracavity field and demon- 
strated that for the optimal duration of the lasing cycle 
(about a few durations of the superradiation pulse of a 
single atomic ensemble), a phase diffusion coefficient D 
as low as a few 10~^ can be realized. We found that this 
phase relay process is robust to 10 % fiuctuations of the 
atom number introduced into the cavity. We conclude 
that the method of sequencial coupling may be a promis- 
ing approach towards the creation of an active optical 
frequency standard. 
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APPENDIX A: SELF-CONTRADICTION OF 
CONVENTIONAL CHOICE OF CORRELATION 
FUNCTIONS 

As it was mentioned in lll CI usually the correspondence 
between operator and c-number variables is based on the 
requirement that products of c-number variables corre- 
sponds to the normally ordered operator variables with 
the normal ordering d"*", M+, Na, Nt, M, a [l^. This 
convention gives the diffusion coefficients 2?^; different 
from (fTOl) - (|23p . Instead, one can obtain (see for 
details): 

22?L = 7a(A/'<?(t)) 

- gr^(t) {M*^{t)A{t) + M^{t)A*{t)) , (48) 

2Vl,j^^2gr^{t){M^{t)A{t)), (49) 

2P^.^. = 2gT^{t) {Mlit)A*{t)) , (50) 

2P^^. =(27afc-7a)(AA,(t)), (51) 

22?^A< = 2I?,^Ai. = 0. (52) 



IV. CONCLUSION 

We proposed to keep the phase coherence between 
pulses of a bad-cavity active optical frequency stan- 



If one tries to use these diffusion coefficients for the 
numerical simulation of Langevin forces, one obtains an 
inconsistency. Indeed, if one introduces the real column 
vector ([55]) . then the covariance matrix D occurs to be 
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1 

At 



( ^aNa - GMM* 









V 







GT 



2 J 4 



r 





li 



A/'a(7afc- y) -Gr 



7a A ^^2-^'+-^*' 



(53) 



after adiabatic elimination of the field variables A, A*. 
This matrix must be positivc-scmidefinite because it is 
a covariance matrix of real random values J- a, Kc[J-m], 
Im[J\A4]- However, this is not true in the general case. 
For example, if "fab = 7a/2, the last diagonal element 
of the matrix U is always negative. The first diagonal 
element is also negative, if G ^ 7a (which corresponds 
to the high-cooperativity regime) and if the atomic po- 
larization is not drastically suppressed by some reasons. 
In any of these cases, the matrix D is not positive- 
semidefinite according to Sylvester's criterion. So the 
conventional choice of correlation functions occurs to be 
self-contradicting. 



APPENDIX B: DURATION AND SHAPE OF 
SUPERRADIANCE IMPULSE 

Consider one single ensemble of Af two-level atoms 
coupled to the resonance mode of the bad cavity in the 
high-cooperativity regime, i.e. when Afg^ S> 7abK. We 
study the time dependence of the intracavity field us- 
ing the set of equations for c- number variables A, Ai, 
Ma when r(<) = 1. For the sake of simplicity we ne- 
glect the Langevin forces and suppose that initially all 
the atoms arc in the upper lasing state. Of course such 
simplification obstructs the correct behaviour when the 
superradiance pulse starts. However, this model can re- 
produce the shape of the pulse because the interaction of 
the atomic ensemble with the thermal bath will play a 
significant role only in the very beginning, or in the very 
end of the pulse. 

We adiabatically express the field variable A = 2gAA / k 
via the atomic polarization AA similarly to (P7)) . Then 
the set of equations for the atomic variables reads as 



Afa 



AA = -"iahAA + GM 



(54) 



-JaMa-GMM*, 



(55) 



where G — Ag^/n. In the high-cooperativity regime we 
neglect the first terms on the right-hand side of these 
equations. Then we introduce Q = A4Ai* . The set of 
equations transforms to 



Q - G{M~2Ma)Q 

Afa = ~GQ. 
The solution to this set of equations reads as 



Q{t) 



Afa{t) 



Af^ 



4cosh^[G7V(<-to)/2] 



1 + cxp[GAf{t - to)]' 



(56) 
(57) 



(58) 
(59) 



where to is the time when the atomic polarization (and 
therefore the intracavity field) is maximal. Here we took 
into account that A/'a(— oo) = Af, A4{~go) = 0. 

The superradiance pulse has the shape 



Ait) 



2g Afexp{i(j)) 
~2cosh[G7V(t-to)/2]' 



(60) 



where (j) is the phase of the intracavity field. The char- 
acteristic duration tp of this pulse (measured to the half 
of its height) is 



4arccosh2 5.27 



GAf 



GN' 



(61) 
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